Impact of hollow-atom formation on coherent x-ray scattering at high intensity 
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X-ray free-electron lasers (FELs) are promising tools for structural determination of macromolecules via 
coherent x-ray scattering. During ultrashort and ultraintense x-ray pulses with an atomic scale wavelength, 
samples are subject to radiation damage and possibly become highly ionized, which may influence the quality 
of x-ray scattering patterns. We develop a toolkit to treat detailed ionization, relaxation, and scattering dynamics 
for an atom within a consistent theoretical framework. The coherent x-ray scattering problem including radiation 
damage is investigated as a function of x-ray FEL parameters such as pulse length, fluence, and photon energy. 
We find that the x-ray scattering intensity saturates at a fluence of ~ 10^ photons/ per pulse, but can be 
maximized by using a pulse duration much shorter than the time scales involved in the relaxation of the inner- 
shell vacancy states created. Under these conditions, both inner-shell electrons in a carbon atom are removed, 
and the resulting hollow atom gives rise to a scattering pattern with little loss of quality for a spatial resolution 
> 1 A. Our numerical results predict that in order to scatter from a carbon atom 0. 1 photons per x-ray pulse, 
within a spatial resolution of 1.7 A, a fluence of 1 x 10^ photons/ A^ per pulse is required at a pulse length of 1 fs 
and a photon energy of 12 keV. By using a pulse length of a few hundred attoseconds, one can suppress even 
secondary ionization processes in extended systems. The present results suggest that high-brightness attosecond 
x-ray FELs would be ideal for single-shot imaging of individual macromolecules. 

PACS numbers: 32.90.+a, 32.80.Fb, 87.59.-e, 87.15.ag 



I. INTRODUCTION 

X-ray free-electron lasers (FELs) [1, 2] provide unparal- 
leled peak brightness and open a new era in science and tech- 
nology, offering many possibilities that have not been con- 
ceivable with conventional light sources [3-5]. The world's 
first x-ray FEL — the Linac Coherent Light Source (LCLS) at 
SLAC National Accelerator Laboratory [6] — has been in op- 
eration since 2009, with a photon energy of up to 8.3 keV, up 
to 2 X lO'^ photons per pulse, and a full-width-half-maximum 
(FWHM) pulse length as short as a few femtoseconds. The 
SPring-8 Compact SASE Source (SCSS) at SPring-8 [7] and 
the European X-ray FEL at DESY [8] are under construction 
and are planned to deliver up to 12 keV photon energy with 
an average brightness 5-500 times higher than LCLS. 

One of the prospective applications of x-ray FELs is single- 
shot imaging of individual macromolecules [9-11], which 
employs coherent x-ray scattering to determine the atomically 
resolved structure of non-crystallized biomolecules or other 
nanoparticles [12-20]. Single-shot imaging becomes possi- 
ble because the high fluence of a tightly focused x-ray FEL 
pulse could produce a significant amount of scattered photons 
from single-molecule samples. One of the key challenges in 
single-shot imaging using ultraintense x rays is radiation dam- 
age [21]. Each target molecule undergoes electronic dam- 
age via processes such as photoionization and Auger decay. 
The positively charged atomic ions formed in this way re- 
pel each other, thus leading to Coulomb explosion of the tar- 
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get molecule [22-24]. Since the fluence required for single- 
shot imaging exceeds the conventional damage limit (200 
photons/A^) [25], these damage effects could degrade the 
scattering patterns and hinder the determination of the atomic 
positions in the target molecule. To suppress the impact of the 
molecular Coulomb explosion on atomically resolved imag- 
ing, one must effectively freeze the atomic motion during the 
x-ray pulse, requiring a pulse duration of no more than ten 
femtoseconds [26, 27]. Following this idea, several theoretical 
approaches were used to simulate the radiation damage pro- 
cesses including the movement of the atoms: References [27- 
33] employed molecular dynamics, Refs. [34-38] based their 
description on a hydrodynamic model, and Refs. [24, 39-42] 
used a kinetic Boltzmann model. 

For electronic damage processes, which in the x-ray regime 
are mainly atom-specific, one may concentrate, to a first ap- 
proximation, on the interaction of the x rays with individual 
atoms. The dynamics of bound electrons in an isolated atom 
during an ultraintense x-ray pulse were investigated theoreti- 
cally in connection with hollow-atom formation [43, 44], x- 
ray fluorescence [45-47], and radiation damage [48]. It is im- 
portant to note that even if the atomic motion during an x-ray 
pulse is negligible, electronic damage dynamics during the x- 
ray pulse may directly influence x-ray scattering patterns by 
altering the electronic density in the target [49]. It is, there- 
fore, crucial to understand detailed ionization and relaxation 
dynamics in individual atoms under ultrashort and ultraintense 
x-ray pulses. 

A series of recent experiments conducted at LCLS re- 
vealed how electrons interact with ultraintense, ultrafast x-ray 
pulses [50-53]. In the x-ray regime, photoabsorption predom- 
inantly affects inner-shell (core) electrons. If inner-shell pho- 
toabsorption is saturated, all inner-shell electrons in a given 
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atom may be removed before Auger decay or other relax- 
ation processes occur [43, 44]. The transient state thus pro- 
duced is referred to as a hollow atom. If the pulse length is 
short enough, the hollow atom retains its core vacancies dur- 
ing the pulse. Because the x-ray photoabsorption probability 
is smaller for valence electrons than for inner-shell electrons, 
hollow-atom formation suppresses further electronic damage. 
This effect is called x-ray transparency [50] or frustrated ab- 
sorption [51], and might be beneficial for single-shot imaging 
of individual molecules [50, 51]. The present paper investi- 
gates this idea in detail and provides criteria for the x-ray FEL 
parameters required. 

To treat x-ray-atom interactions, we employ a consistent 
ab initio framework [54] based on nonrelativistic quantum 
electrodynamics and perturbation theory. This x-ray atomic 
theory has been applied to study x-ray absorption by laser- 
dressed atoms [55-58] and x-ray scattering from laser-aligned 
molecules [59-61]. In this paper, we present a practical imple- 
mentation of this ab initio framework as a toolkit to calculate 
cross sections and rates of x-ray-induced processes for various 
charge states and electronic configurations of an isolated atom 
within an approximation to the Hartree-Fock model. With 
those parameters, we simulate hollow-atom formation dynam- 
ics under ultrashort and ultraintense x-ray pulses by means of 
time-dependent rate equations [44, 50, 62]. Then we investi- 
gate the effects of hollow-atom formation on coherent x-ray 
scattering from atomic carbon. 

The paper is organized as follows. In Sec. II, we present 
theoretical methods to compute cross sections and rates of all 
electronic damage processes, which are integrated into a set 
of rate equations. In Sec. Ill, we explore coherent x-ray scat- 
tering signals influenced by hollow-atom formation and their 
dependence on x-ray FEL parameters such as pulse length, 
fluence, and photon energy, and on the spatial resolution of 
the image. We also discuss the role of electron impact ioniza- 
tion in molecules. We conclude with a summary and future 
perspectives in Sec. IV. 



II. THEORY AND NUMERICAL DETAILS 

The present toolkit of x-ray atomic processes covers ion- 
ization, relaxation (Auger decay and fluorescence), and coher- 
ent x-ray scattering. We assume that inelastically (Compton) 
scattered photons are energetically distinguishable from elas- 
tically (coherently) scattered photons and focus on coherent 
scattering processes for imaging problems. Compton scatter- 
ing contributes to electronic damage, but is negligible in com- 
parison with photoionization for the photon energies under 
consideration [63]. Shake-up and shake-off processes [64- 
66] also make a small contribution to electronic damage and 
are not included in our model. We also neglect impact ion- 
ization [28, 34, 48, 67, 68], i.e., secondary ionization in 
molecules induced by photoelectrons and/or Auger electrons 
via (e,2e) processes. We will discuss a straightforward strat- 
egy to reduce impact ionization in Sec. Ill E. Atomic units are 
used in this section. 



A. Hartree-Fock-Slater model 

In order to implement the ab initio framework [54], we 
use the Hartree-Fock-Slater (HFS) model [69, 70], which 
employs a local density approximation to the exact ex- 
change interaction. The effective one-electron (mean-field) 
Schrodinger equation to be solved is 
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Here the potential is given by 
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where Z is the nuclear charge and the electronic density p (r) 
is given by 
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where / is the spin-orbital index and A^eiec is the number of 
electrons. The exchange term is approximated by the Slater 
exchange potential [69], 
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In addition, the potential includes the Latter tail correc- 
tion [71] to obtain the proper long-range potential for both 
occupied and unoccupied orbitals, i.e., we put V'(r) = — (Z' + 
l)/r if the right-hand side of Eq. (2) is less negative than 
— (Z' + 1 ) /r, where Z' = Z — A^eiec is the effective charge of 
the system. 

After angular momentum averaging, the problem becomes 
spherically symmetric, and each solution of Eq. (1) can be 
expressed in terms of the product of a radial wave function 
and a spherical harmonic. For example, a bound-state spatial 
orbital with quantum numbers {n,l,m) may be written as 



r„/,„(r) = ^y/"(0,</)). 



(5) 



For the bound states, the radial wave function P„i{r) is accu- 
rately solved by the generalized pseudospectral method [72, 
73] on a nonuniform grid. For the continuum states, Peiir) is 
numerically solved by the fourth-order Runge-Kutta method 
for a given energy £ on a uniform grid [74, 75]. To evaluate 
integrals involving both bound and continuum states, we use 
spUne interpolation to map the bound-state orbitals from the 
nonuniform grid to the denser uniform grid employed for the 
continuum states. 

For the bound-state calculations, the theoretical procedure 
presented here is identical to the Herman-Skillman code [70], 
which has been widely used in atomic physics, but the numer- 
ical part of the present toolkit utilizes a different grid method 
with the following advantages. Firstly, it is easy to control 
convergence with respect to the grid parameters. Secondly, 
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we can avoid truncation of the maximum radius internally im- 
posed by the Herman-Skillman code, which causes numeri- 
cal instability for photoabsorption cross section calculations. 
Lastly, the matrix eigenvalue problem is solved by a mod- 
ern linear algebra package [76]. For the present calculations, 
we use 200 grid points and a maximum radius of 50 a.u. for 
the bound states (nonuniform grids) and a radial step size of 
0.001 a.u. for the continuum states (uniform grids) to achieve 
machine accuracy for cross sections and rates. 



2, 3, 4, 5, or all 6 electrons from the occupied orbitals is 
1+3 + 6 + 7 + 6 + 3 + 1 = 27. Note that we perform a separate 
HFS calculation for each configuration. In other words, the or- 
bitals are optimized in the presence of core and/or valence va- 
cancies. Thus, orbital relaxation for the core-hole configura- 
tions is automatically included, a strategy that is known to be 
in good agreement with multiconfigurational self-consistent- 
field calculations [49] . 



We calculate all cross sections and rates for all possi- 
ble electronic configurations. For example, a neutral carbon 
atom has a ls^2s^2p^ ground configuration, so the number 
of all configurations that can be formed by removing 0, 1, 
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B. X-ray absorption process 

The cross section for ionizing an electron in the ith subshell 
by absorbing an x-ray photon with energy (O is given by [54] 
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where a is the fine-structure constant, A^, is the occupation 
number of the ith subshell, /> = max(Z, , Z^), and e = o — £, is 
the photoelectron energy. Here, £, is the ionization energy of 
the ;th subshell (£, = — e,) by Koopmans' theorem [77], which 
is approximately valid in the HFS model. The orbital energy 
£,■ and the radial wave functions P„./. (r) and Pgi. (r) are calcu- 
lated for a given electronic configuration. We do not consider 
orbital hole alignment after ionization by linearly-polarized x- 
ray pulses and, hence, assume that the density of bound elec- 
trons remains spherically symmetric throughout. 

Table I shows x-ray absorption cross sections for all con- 
figurations of carbon (except the bare nucleus) at photon en- 
ergies of 8 keV and 12 keV, respectively. The photon energy 
range considered here is consistent with that available at cur- 
rent and future x-ray FELs [6-8]. Note that the present total 
cross sections are in agreement with those in the literature [78] 
to within less than 10%. 

The impact of orbital relaxation is evident in the subshell 
cross sections shown in Table I. For instance, the 2s subshell 



cross section for the ls'2s^2p^ and ls^2s'^2p^ configurations 
differ by 40% and 90%, respectively, from the 2s subshell 
cross section for the ground configuration of neutral carbon. 

As may be seen in Table I, the cross sections for the Is 
and 2s subshells are much greater than the 2p subshell cross 
section, because Op is proportional to o^'^^/^ in the high en- 
ergy limit [79]. Therefore, absorption of linearly polarized x 
rays does not, in general, induce any orbital hole alignment. 
This justifies our assumption of spherically symmetric elec- 
tron densities. 



C. Auger decay process 

The Auger decay rate that an electron from the yth subshell 
fills the ;th subshell and another electron from the /th subshell 
is ejected into the continuum may be written as [54, 80] 



+ ^ L=\lj-lj,\S=0 /,./ 

where i' indicates the continuum state with Auger electron energy e = £, —Ej — £/, A^,^ is the number of holes in the ith subshell, 
and 

{ /Ai .'I'ul, I o\ for inequivalent electrons , 
(40+2)(4/.,+2) -1 
,.41S;;l„ for =,uiv.le„. deepen. 

Here, averaging schemes over initial and final states to compute transition rates are adopted from Refs. [80-82]. 
The matrix element M^s is given by 

MLsUJ',i,i') - T(-1)^+'^H' £ [RK{jj',i,i')AKUj\i,i') + {-1)'+'RkU' J,i,i')Mj' J,i,i')] (9) 
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TABLE I. X-ray absorption cross sections ((Tp) for various configurations of carbon at 8 keV and 12 keV. 









CTp (10^^ a.u.) at8keV 






(7p (10-*^ a.u.) at 12 keV 




Charge 


Configuration 


is 


2s 


2p 


Is 


2s 


2p 


+0 


ls^2s^2p^ 


287 


14.8 


0.0897 


77.8 


4.05 


0.0164 


+ 1 


ls'2s'^2p^ 
Is22si2p' 


155 


20.8 


0.219 


41.9 


5.70 


0.0399 




287 


8.28 


0.118 


78.0 


2.27 


0.0214 




Is22s22p' 


287 


16.6 


0.0590 


77.9 


4.53 


0.0107 


+2 


ls"2s'^2p- 
Isi2si2p' 


- 


27.9 


0.387 


- 


7.61 


0.0699 




155 


11.5 


0.258 


42.0 


3.16 


0.0474 




Isi2s22p' 
ls22s"2p- 


155 


23.4 


0.132 


42.0 


6.41 


0.0243 




288 




0.145 


78.2 


- 


0.0266 




ls22s'2p' 
Is22s22p0 


288 


9.25 


0.0737 


78.2 


2.53 


0.0134 




288 


18.6 


- 


78.1 


5.10 


- 


+3 


ls"2si2p'= 


- 


15.1 


0.435 


- 


4.13 


0.0790 




ls°2s22pl 
Isl2s02p2 


- 


31.1 


0.221 


- 


8.48 


0.0397 




156 


- 


0.300 


42.1 




0.0548 




Isi2si2pi 
Isi2s22p0 


156 


12.8 


0.153 


42.2 


3.51 


0.0280 




156 


26.0 


- 


42.3 


7.11 


- 




ls-2s"2pi 


289 


- 


0.0885 


78.4 




0.0165 




Is-2si2p0 


289 


10.8 


- 


78.5 


2.93 


- 


+4 


ls"2s"2p^ 


- 


- 


0.454 


- 


- 


0.0841 




Is02si2pi 


— 


16.6 


0.246 


— 


4.54 


0.0442 




1 sO?s29„0 

la jLS jLU 




JJ.O 












Isi2s02pl 


156 




0.180 


42.2 




0.0323 




Isi2si2p0 


157 


14.7 




42.4 


4.03 






ls22s°2p'^ 


284 






76.3 






+5 


ls"2s"2p^ 






0.285 






0.0522 




ls°2s'2p" 




18.0 






4.91 






ls'2s°2p" 


156 






42.2 







where T = 1 / \/2 if j and / are equivalent electrons and 
T = 1 otherwise. Ak is a coefficient related to 3j- and 6j- 
symbols [83], 
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where 
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(/||Q||Z') = (-l)V(2/ + l)(2/' + l)(^0 
and Rk is a double radial integral defined as 
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RKUj',i,i') = Pn,ij{ri)Pnyiy {r2)^Pn,i.{n)Pei,, {r2) dn dr2. 
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Table II lists the Auger rates computed using Eq. (7) for var- 
ious configurations of carbon. KLiL\, KL\L2t,, and KL2t,L2t, 
represent rA(ls,2s2s), rA(ls,2s2p), and rA(ls,2p2p), respec- 
tively. For the Is hole configuration, the present results are 
in good agreement with other theoretical results [47, 81, 84]. 
Note that Ref. [47] includes Auger and fluorescence rates for 
all possible configurations of carbon computed by Cowan's 
atomic structure code [85]. The decay rates in Ref. [47] are 
in fair agreement, to within about 40%, with the present re- 
sults. We note that the experimental Auger lifetime for a free 
carbon ion, averaged over the subset of all doublet states in 
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the ls^2s^2p^ configuration, is 7.3 fs [86], and possibly for- 
tuitously agrees somewhat better with our calculated value 
(wlO fs) than with Ref. [47] («14 fs). 



D. Fluorescence process 

The fluorescence rate for the electric dipole transition of an 
electron from the y'th subshell to a hole in the /th subshell is 
given by [54, 80] 
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The last two columns in Table II show the x-ray fluorescence 
rates and yields [= rF/(rF + ^rA)] for various configura- 
tions of carbon, where Ka indicates rF(ls,2p). For the Is 
hole configuration, the x-ray fluorescence yield is compared 
with experimental data [87]. In light atoms like carbon, the 
x-ray fluorescence yield is generally small. Note, however, 
that x-ray fluorescence is the only decay process available in 
ls'2s*'2p' and ls''2s*'2p^. All fluorescence rates are 
included in our damage dynamics model for completeness. 



E. Rate equations for ionization and relaxation dynamics 

To simulate electronic damage dynamics in intense x-ray 
pulses, we use the rate equation approach that has been 
successfully used to describe x-ray-induced multiple ioniza- 
tion [44, 50, 62]. The transitions among all possible electronic 
configurations {/} of a given atom are described by a set of 
coupled rate equations of the form 



all config 



at 



(14) 



where Pj is the population of the /th configuration, and F/^// 
is the rate for transitions from the configuration / to the con- 
figuration /'. Here F can be either a time-independent Auger 
or fluorescence rate, or a time-dependent photoionization rate 
given by apJ{t), where 7(f) is the photon flux of the x-ray 
pulse at a given time t. In our calculations on carbon, all 
configurations connected by the photoionization. Auger de- 
cay, and x-ray fluorescence processes listed in Tables I and II 
are included, and the corresponding rate equations are numer- 
ically solved using the fourth-order Runge-Kutta method. We 
assume that the temporal shape of the x-ray pulse is Gaussian. 
In the regime considered here, the spiky structure of the in- 
dividual pulses generated by x-ray FELs such as the LCLS is 
largely irrelevant [44] . 



F. Coherent x-ray scattering process 

The coherent x-ray scattering form factor for a given elec- 
tronic density p (r) is given by [54] 



fiQ)= /p(r).'Q-'^V, 



(15) 



where Q is the photon momentum transfer We assume that 
the atomic electron density is spherically symmetric. Then 
the atomic form factor depends only on the magnitude of the 
momentum transfer, so Eq. (15) may be simplified to 



Jo Qr 
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where 2 = |Q| = 2aa)sin(0/2) and Q is the polar angle of the 
momentum of the scattered photon with respect to the propa- 
gation axis of the incoming x rays. 

For the unpolarized x rays, the differential cross section for 
coherent scattering is given by 



(17) 



and for hnearly polarized x rays, the differential cross section 
is given by 



«^I/"(G)|'(1 -cos^^sin^e), (18) 



where is the azimuthal angle of the scattered photon mo- 
mentum with respect to the x-ray propagation and polariza- 
tion axes. This differential cross section gives the x-ray scat- 
tering pattern one would obtain for a fixed electronic configu- 
ration. From measurement of the x-ray scattering pattern, one 
can retrieve electronic density information. Based on the dis- 
crete Fourier transform relationship between real space and 
Q space, the smallest length in real space (the resolution d) 
corresponds to the largest length in Q space (Qmax), while the 
largest length in real space (the object size D) corresponds to 
the smallest length in Q space (the pixel size AQ). In other 
words, the spatial resolution desired determines the photon 
momentum transfer up to which statistically significant scat- 
tering data must be available, Qmax = 2n/d; the size of the 
object in real space determines the maximum pixel size per- 
mitted in Q space, Ag = 2k/D. For a purely atomic target 
(the case considered here), the object size is close to the desir- 
able spatial resolution, so all relevant information in momen- 
tum space is captured in one pixel. The number of photons 
scattered into that pixel is proportional to the integral of the 
differential cross section over the solid angle i2 up to the de- 
sired resolution. We note that the integrals of Eqs. (17) and 
(18) over H are identical. 

During exposure to an ultraintense x-ray pulse, the atomic 
electron density is dynamically modified, as a consequence 
of x-ray-induced processes. This makes it necessary to intro- 
duce a suitably averaged, time-dependent differential scatter- 
ing cross section. 



— it) 
dO. ^ ' 



all config. 



(19) 



(16) 



where P/(f) is the population of the /th configuration, which 
is obtained as a solution of the rate equations in Eq. (14). The 
differential scattering cross section for the /th configuration 
is evaluated from the form factor of Eq. (16) using the den- 
sity p(r) calculated from the orbitals optimized for the /th 
configuration. In this way, we incorporate electronic damage 
dynamics into simulations of coherent x-ray scattering at high 
intensity. 
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TABLE II. Auger rates (Fa) and fluorescence rates (Fp) for various configurations of carbon. For the Is hole configuration, the present 
values are compared with other approaches. A: Herman-Skillman code [84], B: Cowan code [47], C; semi-empirical method [81], and EXP: 
experimental data [87]. 









Fa (10-3 a.u.) 




Fp (10-5 ^ y ) 


Fluorescence 


Charge 


Configuration 


KLiLi 






Ka 


yield 


+1 


ls'2s'^2p^ 














Present 


0.961 


0.970 


0.439 


0.836 


0.0035 




A 


0.929 


0.987 


0.435 


0.824 


0.0035 




B 


0.680 


0.697 


0.392 


0.651 


0.0037 




C 


0.857 


0.824 


0.378 


0.486 


0.0024 




EXP 


— 


— 


— 




0.0026 


+2 


ls"2s'^2p^ 


2.89 


3.33 


1.75 


■ : 


0.0034 




Isi2si2p2 


— 


0.602 


0.574 


0.975 


0.0082 




Isi2s22pl 


1.18 


0.620 


— 


0.498 


0.0028 


+3 


Z — /\_ 1 _ — T 

ls''2si2p' 




1.98 


2.13 


3.04 


0.0073 




ls°2s22pi 


3.46 


1.99 




1.55 


0.0028 




Isi2s°2p2 






0.703 


l.ll 


0.0155 




Isi2si2pl 




0.370 




0.569 


0.0151 




Isi2s22p° 


1.39 










+4 


ls"2s"2p- 






2.59 


3.27 


0.0125 




Is02si2pl 




1. 15 




1.69 


0.0145 




ls°2s22p'' 


3.91 












lsi2s''2p' 








0.648 


1.0000 


+5 


ls"2s''2p' 








1.97 


1.0000 



Measurement of the strength and the quality of scattering 
signals 



To describe the strength of the x-ray scattering signals, we 
compute the number of scattered photons (A^s) by integrating 
over the Ewald sphere [88], limited to the desired spatial res- 



olution length d, 

Ns{d) = 



e<eraax(f/) 



(20) 



where 7(f) is the incident photon flux at a given time t, and the 
time-dependent differential scattering cross section is defined 
in Eq. (19). Here Qmax is determined by the spatial resolution 
length, d = 27i;/gmax- 



To measure the quality of the x-ray scattering patterns, we employ a modified /^-factor expression with an explicit dependence 
on the spatial resolution [35], 



Rid) 



v/jVreal(^^) 



VMdeal(^^) 



,/ v/Mdeal(ii') dQ.' 

e'<eraax(^/) 



da. 



(21) 



where 

is the number of photons (per unit solid angle) scattered from 
the sample undergoing x-ray-induced electronic damage, and 

Nii,.xm^( r J{t)dt\^ =F^ (23) 

\J-oo / CtLl neutral dii neutral 

is the number of photons (per unit solid angle) scattered from 
the undamaged sample, which is given by the fluence (F) 



times the differential cross section of the undamaged sample. 
In our case, the undamaged sample is the neutral carbon atom 
in its ground configuration. 



III. RESULTS AND DISCUSSION 

We investigate high-intensity coherent x-ray scattering in- 
cluding electronic damage dynamics with x-ray parameters 
achievable using x-ray FELs [6-8]. The FWHM pulse length 
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FIG. 1 . (Color online) Normalized atomic form factor for different 
electronic configurations. 



in our calculations varies from 1 fs to 120 fs, the pulse 
envelope being a Gaussian. The number of incident pho- 
tons varies from 10^ to lO'^, and the beam size used is 
100 nmx 100 nm, corresponding to a fluence ranging from 
10^ to lO'' photons/A^. The photon energy is chosen as 8 keV 
and 12 keV, respectively. With the Gaussian envelope and 
these fluences, the peak intensity ranges from 1 x lO'^ to 
1 X 10^^ W/cm^, for a photon energy of 8 keV and a pulse 
length of 120 fs. 



A. Atomic form factors for core hole states 

In order to examine variations of x-ray scattering patterns 
for different electronic configurations, especially for core- 
hole configurations created via photoabsorption, we calculate 
atomic form factors for the filled core (neutral: ls^2s^2p^), 
the single-core-hole (ls'2s^2p^), and the double-core-hole 
(ls''2s^2p^) configurations (see Fig. 1). To facilitate a direct 
comparison between the three different charge states, the form 
factors in Fig. 1 are normalized in the same fashion as used in 
the /^-factor expression in Sec. II G. Specifically, the normal- 
ization factor is 



e<eraax(^/) ' 



^NiiQ.) d£l, (24) 



e<emax(f/) 



where A^/(i2) = das/dQ. for a given configuration/. We keep 
J = 1.7 A fixed as a desirable resolution for further analy- 
sis. Within this resolution, i.e., for Q up to Qmax, the shapes 
of the three normalized form factors are quite similar to each 
other. The computed /^-factors for the single-core-hole and 
double-core-hole configurations with respect to the ground 
configuration of the neutral atom are 1.7% and 2.6%, respec- 
tively. This fact indicates that core-hole formation causes little 
degradation of the quality of the x-ray scattering pattern. We 
will discuss the optimal resolution minimizing the /J-factor in 
Sec. HID. 



B. Hollow-atom formation in ultrashort and ultraintense x-ray 

pulses 

The time-averaged charge weighted by the normalized 
pulse envelope provides a simple measure of electronic dam- 
age during the x-ray pulse [48, 49]. The time-averaged popu- 
lation of the /th configuration is given by 



(25) 



Pi^ / Pi{t)f{t)dt, 



where Pi{t) is the time-dependent population of the /th con- 
figuration, and /(f) is the normalized x-ray pulse envelope. 
Then the time-averaged charge is given by 

all config. 

z= Z'P'^ 
I 

where Z/ is the charge corresponding to the /th configuration. 

Figure 2 shows the time-averaged charge of atomic carbon 
as a function of the pulse length for several fluences. When the 
pulse length is short enough to compete with core-hole life- 
times as marked in Fig. 2 (Ts « 10 fs for the single-core-hole 
configuration and Td ~ 3 fs for the double-core-hole configu- 
ration), Z starts to decrease. This reduction of electronic dam- 
age is the signature of x-ray transparency [50] or frustrated 
absorption [51], which may be understood as follows. Pho- 
toionization of a core electron initiates electronic damage. If 
the x-ray pulse length is long enough for Auger decay to occur 
during the pulse, then a valence electron fills the core vacancy. 
Eventually, many electrons can be stripped off in a sequence 
of core photoionization and Auger decay steps. Note that in 
this limit, the effective x-ray absorption cross section remains 
essentially constant throughout the pulse. For instance, the x- 
ray absorption cross sections for the configurations ls^2s"2p^ 
(neutral ground configuration) and ls"2s^2p*' (the dominant 
Auger decay channel for the single-hole configuration) are al- 
most identical. As Fig. 2 illustrates, the time-averaged charge 
is practically independent of the pulse duration (and, there- 
fore, independent of the peak intensity) for pulse lengths much 
longer than 10 fs. For pulses shorter than 10 fs, there are two 
stages of x-ray transparency or frustrated absorption. In the 
first stage, if the pulses are still longer than the double-core- 
hole lifetime (Td ~ 3 fs) and only single-core-hole production 
is saturated, the effective x-ray absorption cross section drops 
by a factor of 1 .7 relative to the neutral atom (cf. Table I). The 
second stage becomes accessible for pulse durations shorter 
than the double-core-hole lifetime. In this case, by saturating 
hollow-atom formation, the effective x-ray absorption cross 
section drops further by a factor of 6.2, i.e., relative to the 
neutral atom the effective x-ray absorption cross section drops 
by almost a factor of 1 1 (cf. Table I). The somewhat counter- 
intuitive consequence of this is that by decreasing the pulse 
duration and, thereby, increasing the peak intensity, the time- 
averaged charge Z can be minimized, as shown in Fig. 2. 

We compare the time-averaged charge for two different 
photon energies in Fig. 2: (a) 8 keV and (b) 12 keV. For a 
given fluence, the time-averaged charge at 8 keV is higher 
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FIG. 2. (Color online) Time-averaged charge weighted by the normalized pulse envelope for 8 keV and 12 keV photon energies. Each line 
corresponds to a different fluence F (photons/A^). The vertical lines labeled by Ts and Td indicate the single-core-hole and double-core-hole 
lifetimes, respectively. 



than that at 12 keV, because the photoabsorption cross sec- 
tion at 8 keV is about 4 times higher than that at 12 keV (Ta- 
ble I). This is also expected based on the scaling behavior 
of the photoabsorption cross section in the high energy limit, 
C7p o<: O^^/^ {I = 0) [79]. Since a higher photon energy in- 
duces less electronic damage, it has an advantage with respect 
to the /^-factor for x-ray scattering, which will be discussed in 
Sec. HID. 

As one can see in Fig. 2, Z increases in increments of de- 
creasing magnitude as the fluence becomes higher, indicat- 
ing a saturation effect. To make this point clearer, we plot 
the time-averaged population P/ of the single (/=ls'2s^2p^) 
and double (/=ls''2s^2p^) core-hole configurations as a func- 
tion of the fluence in Fig. 3. The pulse length is fixed at 1 fs 
FWHM, which is less than the lifetimes of the core-hole con- 
figurations. Both populations follow a power-law dependence 
for fluences up to about 10'' photons/ A^: single-core-hole pro- 
duction is a one-photon process and therefore is a linear func- 
tion of the fluence (below 10^ photons/ A^), whereas double- 
core-hole production is a two-photon process and therefore is 
a quadratic function of the fluence. Saturation occurs around 
10^-10** photons/A^. For even higher fluences, the popula- 
tions of both configurations decrease due to further photoion- 
ization (core ionization in the single-core-hole configuration 
and valence ionization in the double-core-hole configuration, 
respectively). 



C. Influence of hollow-atom formation on x-ray scattering 
intensity 

The pulse-length and fluence dependence of the time- 
averaged charge affects the number of scattered photons, 
which must be maximized in single-shot experiments such as 
to obtain an optimal signal-to-noise ratio. In Fig. 4, the num- 
ber of scattered photons is plotted as a function of the fluence 
for 8 keV and 12 keV photon energies. The spatial resolution 



is fixed at ii= 1 .7 A and three different pulse lengths are consid- 
ered(l, 10, and 100 fs). For fluences below ^ 10^ photons/A^, 
A^s depends linearly on the fluence of incident photons, but 
is independent of the pulse length (see Fig. 4). In this low- 
fluence regime, the number of photons scattered per atom and 
per pulse is less than 0.1. In order to scatter at least 0. 1 pho- 
tons, the fluence must be in the regime above 10^ photons/A^. 
At high fluence, after saturation of inner-shell ionization, A^s 
is no longer a linear function of the fluence and, particularly, 
depends sensitively on the pulse length. As can be seen in 
Fig. 4, the number of photons scattered may be maximized 
at a given fluence by using a pulse duration shorter than the 
double-core-hole lifetime (Td ~ 3 fs). 

If we require that a carbon atom scatters, say, 0.1 pho- 
tons per pulse and per pixel, then, assuming a pulse length 
of 1 fs and a photon energy of 12 keV, a fluence of 1 x 10^ 
photons/A^ per pulse is needed. However, if we assume a 10 
fs pulse instead, then the fluence required would increase by 
a factor of four. Figure 4 illustrates quite distinctly the impact 
of hollow-atom formation on coherent x-ray scattering at high 
intensity. In a molecule consisting of A^atom atoms, the num- 

1/3 

ber of scattered photons is proportional to at least N^l^^ per 
pulse and per pixel [89]. For example, with the above x-ray 
parameters, a molecule consisting of 100,000 carbon atoms 
would scatter at least 5 photons per pulse and per pixel. Note 
that 5 photons per pixel would be sufficient for successful 3- 
dimensional structural reconstruction [90]. 



D. Dependence of the i?-factor on the desired resolution 

In addition to the strength of the scattering signal, another 
important factor is the quality of the x-ray scattering pattern. 
The scattering pattern from the damaged sample should be 
as similar as possible to the scattering pattern that would be 
obtained if the sample were unaffected by radiation damage. 
Using the /^-factor in Eq. (21), we measure the quality of the 
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FIG. 3. (Color online) Time-averaged populations of the single-core-hole and double-core-hole configurations for 8 IceV and 12 keV photon 
energies. The pulse length is fixed at 1 fs FWHM. 
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FIG. 4. (Color online) Number of scattered photons for 8 keV and 12 keV photon energies. A spatial resolution of 1.7 A is assumed. The thin 
dotted lines extrapolate the linear fluence dependence valid at low fluence. 



x-ray scattering pattern. 

In Fig. 5, we examine the correlation between the /^-factor 
and the desired spatial resolution for 8 keV and 12 keV pho- 
ton energies. The pulse length assumed is 1 fs FWHM and 
the fluence is fixed at 10^ photons/A^. Under these condi- 
tions, inner-shell ionization is saturated (see Fig. 3), and ev- 
ery ten pulses about one photon is scattered per carbon atom 
(see Fig. 4). Each curve in Fig. 5 ends at the finest spatial 
resolution possible at the respective photon energy. Because 
of the reduction of electronic damage at higher photon ener- 
gies, the 7J-factor at 12 keV is less than that at 8 keV. For a 
spatial resolution li > 1 A, the spatially localized reduction 
of electron density in the 1 s shell in the single-core-hole and 
double-core-hole configurations is difficult to resolve, render- 
ing the /^-factor rather low. The small local maximum around 
2.5 A is due to the 2s vacancy formed by valence ionization in 
the double-core-hole configuration. For very fine resolution 
{d <\ A), the core vacancy can be resolved, so the /^-factor 




^ Resolution (A) ^° 

FIG. 5. (Color online) iJ-factor as a function of the spatial resolution 
for 8 keV and 12 keV photon energies. The fluence is 10^ photons/ 
and the pulse length is I fs FWHM. 
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rapidly increases. Note, however, that the /^-factor values in 
Fig. 5 are still less than 20%, which is a typical value for x- 
ray crystallographic data [35]. For d= 1.7-1. 9 A, the 7?-factor 
takes on a local minimum value of less than 2%. 



E. Role of impact ionization 

So far, we have focused on electronic damage processes 
that are not strongly affected by the molecular environment. 
An important damage mechanism characteristic of extended 
molecular systems is impact ionization by (quasi-)free elec- 
trons [28, 34, 48, 67, 68]. For an x-ray pulse much shorter 
than inner-shell decay lifetimes, impact ionization by Auger 
electrons is irrelevant for the formation of electronic damage 
during the pulse. On the other hand, impact ionization by 
photoelectrons is not, in general, negligible. Here we discuss 
how to reduce photoelectron impact ionization by using short 
pulses. 

The mean free path for a 12 keV photoelectron in a carbon- 
based medium (diamond) is about 13 nm [68]. Let us assume 
that the photoelectron travels in the x direction in a homoge- 
neous sample. From the definition of the mean free path, it 
follows that the impact ionization probability is given by 

hmMx)^\-e-'l^, (26) 

where x is the distance traveled and A is the mean free path. 
If we allow an impact ionization probability of 20%, the max- 
imum X permitted is — Alog(l — 0.2) = 22% x A. Therefore, 
for molecules with a diameter of 3 nm or less, the impact ion- 
ization probability per photoelectron is less than 20%. For 
much larger molecules, the role of impact ionization can be 
reduced by using an x-ray pulse that is so short that it is over 
before impact ionization has taken place with a probability 
less than a certain percentage. In analogy to Eq. (25), one can 
define the time-averaged impact ionization probability during 
the x-ray pulse as 

^impact = ^ Pimv^a{x{t))f{t)dt, (27) 

where T is the pulse duration and /(f) is the normalized x- 
ray pulse envelope. It is assumed that the electron starts to 
travel at the beginning of the pulse, i.e., x{t) = vt, where 
V is the speed of the photoelectron. If we use a flat-top 
pulse envelope, f{t) = 1/t for < f < T, then Pimpact = 
1 + A/(vT)[exp(— vt/A) — 1]. Figure 6 plots T versus Pimpact 
for a mean free path of 13 nm and a photoelectron energy 
of 12 keV. For a pulse-weighted impact ionization probability 
of 20%, the pulse duration required is about 100 attoseconds, 
corresponding to a Fourier-limited bandwidth of the order of 
10 eV. In calculations using a Gaussian pulse envelope, we 
obtained very similar results. 

IV. CONCLUSION 

In this paper, we have investigated electronic damage and 
coherent x-ray scattering using ultrashort and ultraintense x- 
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FIG. 6. Plot of the pulse duration required for a given impact ioniza- 
tion probability, for a photoelectron with a kinetic energy of 12 keV 
in a carbon-based medium with a mean free path of 13 nm. 

ray pulses attainable with current and future x-ray FELs. For 
all possible electronic configurations of the atomic system, 
we have calculated rate parameters for x-ray-induced damage 
processes including photoionization. Auger decay, and fluo- 
rescence, in a consistent ab initio framework. The impact 
of electronic damage has been studied by employing a set 
of coupled rate equations, which we have incorporated into 
simulations of coherent x-ray scattering signals. We have im- 
plemented an integrated toolkit, XATOM, to treat all above- 
mentioned processes. 

Our numerical simulations of coherent x-ray scattering sig- 
nals including electronic damage dynamics show that hollow- 
atom formation and the associated phenomenon of x-ray 
transparency or frustrated absorption play a crucial role in 
optimizing the strength and quality of single-shot x-ray scat- 
tering signals. Hollow-atom formation is particularly impor- 
tant when the x-ray pulse length is a few femtoseconds or 
shorter, and saturates, in the case of carbon, around a fluence 
of 10^-10^ photons/ A^, corresponding to lO'^-lO'"* photons 
per pulse at a beam size of 100 nmx 100 nm. At a fluence of 
10^ photons/A^, for instance, the number of photons scattered 
per pulse, within a spatial resolution of 1.7 A is about 0.1 per 
carbon atom, at a pulse length of 1 fs and a photon energy of 
12 keV. A hollow atom is resistant to further electronic dam- 
age via photoionization and, for a spatial resolution d> I K, 
gives rise to an x-ray scattering pattern that differs little from 
that obtained for the neutral atom in its ground electronic con- 
figuration. A comparison between our data for 8 keV and 
12 keV photon energies shows that there are no qualitative 
differences. By using a higher photon energy, the quality of 
the scattering pattern, as defined by the P-factor, can be in- 
creased (by reducing electronic damage), but the number of 
photons scattered per pulse decreases somewhat. 

Finally, we have analyzed the role of impact ionization in 
molecules and provided a simple estimate of the pulse dura- 
tion required to suppress impact ionization during the x-ray 
pulse. This estimate, in combination with the calculations pre- 
sented in this paper, suggests that attosecond x-ray FELs [91- 
94] with a pulse length of 100 as, ~' lO'^ photons per pulse. 
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and a photon energy of '-^ 12 keV are ideal for single-shot 
imaging of individual macromolecules at atomic resolution. 
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